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Abstract 

We consider a RG approach for the plasma of magnetic monopoles of the Ioffe- 
Larkin approach to the t-J model. We first derive the interaction parameters of 
the 2+1 plasma of magnetic monopoles. The total charge along the time axis is 
constrained to be zero for each lattice plaquette. Under the one-plaquette approx- 
imation, the problem is equivalent to a one dimensional neutral plasma interacting 
via a potential V(t) ~ t a , with a = 1/3. The plasma is in a dipolar phase if a > 1 
and a possibility of transition towards a Debye screening phase arises if a < 1, 
so that there exists a critical Fermi wave vector k*j such as the plasma is Debye 
screening if kf < k*j- and confined if kf > k*j. The 2+1 dimensional problem is 
treated numerically. We show that k*^ decreases and goes to zero as the number 
of colors increases. This suggests that the assumption of spin-charge decoupling 
within the slave-boson scheme is self-consistent at large enough values of N and 
small enough doping. Elsewhere, a confining force between spinons and antiholons 
appears, suggesting a transition to a Fermi liquid state. 
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1 Introduction 

The physics of strongly interacting electron systems has received considerable attention over 
the recent years, and it still bears many challenges, especially on the theoretical side. Among 
the various methods and ideas which have been explored in this context, gauge theories seemed 
to offer a rather attractive approach 0-0. Their essence is to focuss on the presence of a non 
double occupancy constraint, which leads to a local U(l) symmetry if a slave boson represen- 
tation is used 0-0]. Although it is possible to derive these gauge theories from an expansion 
around a large N mean field theory of the t-J model j4| j5|, they could also be regarded as 
promising candidates for an effective low-energy theory in order to describe for instance the 
anomalous normal state properties of high-T c superconductors. They seem to predict a phase 
diagram for the single band t-J model which qualitatively ressembles the experimental ones 
for copper-oxide superconductors |J. Further more, they reconcile the existence of a large 
Fermi surface corresponding to Luttinger's theorem, as shown by photoemission experiments, 
and the anomalous transport properties, which are mostly governed by holes || 0. Thermo- 
dynamic properties have also been investigated, and a good agreement with high temperature 
expansions for the t-J model has been reached || . However, this work has also pointed out that 
fluctuations of the gauge field are large, in the sense that the variance of the local statistical 
flux around a given plaquette is not small in units of 2n, even down to low temperatures. This 
feature suggests that the presence of the lattice may not be inessential, since it induces a peri- 
odic action as a function of the time and space dependent flux per plaquette. As demonstrated 
by Poliakov, this periodic nature of the gauge field has dramatic consequences on 2+1 dimen- 
sional electrodynamics since it allows for non trivial space time configurations of the gauge field 
(monopoles), which induce charge confinement 0. It should be emphasized that in the context 
of the t-J model, the gauge field Lagrangian density is not the usual —\F iW F iiV term, but it is 
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generated upon integrating out fermionic and bosonic fluctuations || . A perturbative estimate 
of a single monopole action has also been derived in ||, and was found to diverge. However, 
it is clear that some globally neutral configurations (i.e. with the same number of instantons 
and anti-instantons) have a finite action, and the main question is whether the corresponding 
two-component plasma exhibits Debye-screening or not. This viewpoint has been developed 
by Nagaosa |IJ where he assumed a dissipative-type action for the gauge field, which may be 
relevant for the t-J model at small temperatures, since it requires a finite dc conductivity for 
the fermions. His main conclusion is that no major instanton effect is present in the t-J model 
range of parameter since the dissipative nature of the gauge field dynamics strongly inhibits 
quantum tunneling. In the present paper, we address this question from a slightly different 
perspective, with emphasis on the possible zero temperature transitions. By contrast to the 
results of reference [pi, we find that assuming a Ioffe-Larkin form for the monopole plasma 
action leads to a phase transition between a Debye screening phase (which corresponds to a 
confining force between spinons and antiholons), and a dipole phase (leading to unconfined 
spinons and antiholons). The control parameters are the band filling of the underlying t-J 
model, and the number N of fermion colors (the physical case being N = 2). In rather good 
agreement with physical intuition, the dipole phase of the plasma is found at large N and 
small doping. In this regime, spin-charge separation may then be a self-consistent hypothesis. 
We note however that this leads either to a renormalized Fermi liquid or an anomalous liquid 
depending on whether Bose condensation of ho Ions occurs or not. In the other phase, the 
gauge field cannot be treated perturbatively, and the corresponding mesons (bound states of 
spinons and antiholons) are physical electrons. By contrast, a transition to the dipole phase is 
obtained in reference |10| in the presence of an infinitesimal dissipation. The main difficulty of 
the problem is the determination of the phase diagram of the monopole plasma which exhibits 
long range interactions (in space and imaginary time). Furthermore, unlike the case of the 
standard compact 2+1 electrodynamics, a strong anisotropy exists between space and time 
directions, reflecting the lack of Lorentz invariance in the model. This reduced symmetry 
increases the difficulty of real space RG analysis since the functional form of the monopole 
interaction potential is not stable under a RG transformation. Our approach has attempted 
to take advantage of the fact that the interaction is much stronger along the time direction, 
with a r 1//3 dependence. The corresponding one dimensional problem exhibits a phase transi- 
tion between a Debye-screening phase and a dipole phase. We argue and give some numerical 
indication that the unbinding of the monopole-antimonopole pairs along the time direction 
triggers a 2+1 dimensional unbinding, leading to a globally Debye-screening phase. The paper 
is organized as follows: section § defines the statistical problem of the monopole plasma. The 
next section focusses on the 0+1 dimensional problem along a time direction, giving strong 
arguments in favor of a phase transition. The extension to the 2+1 dimensional situation is 
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then discussed, leading to a phase diagram as a function of N and fermion filling, which is 
the main result of the paper. The conclusion is dedicated to a comparison with previous work 
and stresses open questions. 



2 Statistical mechanics of the monopole plasma 

As already stressed in the introduction, we shall assume a Lagrangian of the form 

(r)^c iCT (r) + b la (r)^-b ia (r)) - t, £ (e^&Cj. + h.c.) (1) 




or or 
'b+bj + h.c.) +iJ2 x i (cJ c i- + h Pi - f ) ■ 



The fields Cut(t) and 6j(r) are respectively fermionic and bosonic, and they are defined on a 
two dimensional square lattice with continuous imaginary time. The hopping constants tf and 
tb can be derived from a large N saddle point approximation of the one band t-J model || [§]. 
We shall from now on focus on the effective dynamics of the U(l) gauge field (a^, Aj), assuming 
that fermions and bosons have been traced out. As shown in the references || @, the gauge 
field action to gaussian order is dominated by the fermion contribution at low doping, and 
with the assumptions that the holons have not condensed. Keeping only the transverse part 
which is responsible of the non-Fermi liquid behavior gives || 

( 1 / \ / kk- \ 
S e ff(a,X)=T / - Ui(k,uj)uj 2 + fi(k,u)k 2 ) ld i:j % ^-ai{k,uo)aj{-k, -u) . (2) 



In this equation, 



for \uj\ 2tfkfk, k <C kf and /j,(k,u) = tf/Ylir. The gauge field variables a r<r+n , where n is a 
lattice vector and r a lattice site are denoted in the continuum limit by a n (r + n/2), in order 
to define the two component field aj(p). The time component of a is identified with A. Most 
of the time, we shall use the axial gauge ao = 0. Latin indices such as % and j denote spatial 
components, whereas Greek indices correspond to arbitrary components. The quantities e\ 
and \x are derived with the approximation of a circular Fermi surface, and by taking the long 
wavelength, small frequency limit of the fermion current-current correlation function. The 
main drawback of this action is that the fundamental periodicity of the original action ([]]), 
namely its invariance under — > + 27T is lost. This periodicity allows for non-trivial space- 
time configurations of the field corresponding to tunneling events where the flux threading a 
given plaquette may change by integer multiples of 2ti. Ioffe and Larkin suggest to express (0) 
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in terms of gauge invariant field strength F^ v = d^a u — 9„a M , with /x = 0, 1, 2, and to replace 
the flux par plaquette F 12 by its value modulo 2ir. This algorithm sounds quite natural on 
physical grounds. However, (|]) has been derived perturbatively for 'flat' configurations which 
satisfy Faraday's law: <9 M 6 M = 0, where 

b[i = —£fj,v P F up . (4) 
After the field strength 6 M is taken modulo 27r, it satisfies 

i 

where are the integer charges located at ri,T{. An ambiguity arises in extending the result 
(H) to non trivial configurations. We may add to equation (fj) any quadratic form 

T E / j0j* C ( k > w ) (u 2 b(k,uj)b(-k,-Lu) - k 2 e ± (k,u)e ± (-k,-uj)) (6) 

without changing the result on 'flat' configurations, but the action for non trivial configurations 
will depend on the kernel C(k,u). In @, e±_ and b denote the transverse part of the electric 
field, and the magnetic field respectively. We also note that the perturbative evaluation of the 
fermion loop generates only the function e\ (k, oj)oj 2 +ji{k, u)k 2 . Physical intuition suggests that 
n(k, u>) is identical to the static diamagnetic susceptibility in the u — > limit. This determines 
the two functions ei(k,u) and n(k,u) as given above, and with such a determination, (fj) 
becomes 

r d 2 k 1 1 
S eff (a } \) = TJ2 J j^^^i(k,uj)e ± (k } uj)e ± (-k, -u) + -fi(k,u)b(k,u)b(-k, -u). (7) 

Equation tfjty is then extended to non-trivial configurations thus lifting the ambiguity in the 
choice of the kernel C(k,u). But if the procedure seems perfectly sound at low frequencies, 
the separation between the electric and magnetic parts is less obvious to access at higher 
frequencies. In the bulk of this paper, we assume that this procedure is valid. The action for 
a many monopole configuration with a topological charge q(r, r) = ±2tt is then given by 

s plasma = -E J T^^^ ^^^ + ^^p ^-^)- (8) 

where q(k,u) is the Fourier transform of the charge density, namely 

?(*,«)= f drY,e~ i{k - r+ " T) q(r,T), (9) 

J v> 
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where r is a lattice site. More specifically, this gives 

c Nt frX^ ( d2k nm 

plasma ^ W M (M + £*3) " liUj 

In this formula, the global factor A" has been added. It is simply the number of fermion 
colors in the large TV approaches. Rescaling energies and frequencies by setting tf = ±1, the 
model depends only on two dimensionaless parameters, A" and kf. Assuming a circular Fermi 
surface, the maximal value of kf corresponds to 1/2 electron per site for a given color, which 
gives k f < (2ti) 1 ' 2 . 

Before going further, we should mention that equation (|8|) is not the only candidate for 
the monopole plasma action. Developing an analogy with Josephson junction arrays, and 
emphasizing the dissipative nature of the gauge field, Nagaosa has also considered the following 
action |I0| : 

S d iss = ^- dr dr}^- (1 - cos {a x (r, n, r) - a±(r, n, r ))) . (11) 

47T J—oo JO r n yT — T ) 

for the dissipative part of the gauge field dynamics. As shown in reference [|ll[] , it is possible to 
map it into a statistical model in some regimes, but mostly a bidimensional model is obtained. 
The key variables are the winding number m(r, n) along the time direction: 

a±(r, n, (3) — a±(r, n, 0) = 2nm(r, n) + z/(r, n), (12) 

with z/(r, n) g] — 7T, 7r] and m(r, n) and integer. It seems that both approaches respect the 2ir 
periodicity and the quadratic expansion of the gauge field around a = 0. In the absence of a 
fully first principle derivation, we shall adopt equation (^) as a working hypothesis, and hope 
to clarify this issue in a future work. 

Going back to equation (|ToD , it is important to stress that for any k value, the u integral 
diverges if lim u ^ Q q(k,u) is non vanishing. Therefore, we shall impose a constraint on the 
allowed topological charge configurations, namely that q(k, uj = 0) = for any k. In real 
space, it means that 

/ drq{r, r) = (13 
J o 

for any plaquette located at r. The partition function of the plasma is then 



+ °° 1 2n / fl 3 dr- \ ( 1 

Z -£W5 U ^Sjx(n,...,r„)exp kE«V(. 



r^Ti-Tj)). (14) 



In this equation, r^, Tj denote the space-time coordinates of the monopoles with topological 
charge g^. We set g« = 2tt if 1 < i < n and qi = —2it if n + 1 < i < 2n. X2n( r i, ■■■, r 2n) 
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expresses the constraint and x = 1 if for any r we have 



2n 



J2l^r,ri =0. (15) 



i=l 



The interaction potential v(r,r) is obtained by Fourier transform of equation ( |10D 

, , N „ I d 2 k e" tJ+ " r) . . 

T) = I2^ T 5 1 WfWM^WY (16) 



with 7 = l/6fc/ (we set t/ = 1). An important ingredient in (|14]) is the imaginary time scale To 



which is obtained by calculating the ratio of the two gaussian determinants in the presence and 
in the absence of an instanton. We have carried out this calculation for the broken parabola 
model of Ioffe and Larkin, which leads to 

- sfcn [ tT.I , ? ' . - ) 1 " 

TO 



1 / f d 2 k k 2 \ 

- = ^[tj:j {2n)iei{ktU)U 2 +flk 2) • ( i7 ) 



The interested reader will find a derivation of this result in the appendix. We note that the 
integral is divergent at large frequencies. This may be another signal that we have not yet 
found a satisfactory derivation of this monopole plasma action. Since this To depends on the 
full non linear action of the gauge field, which is still unknown, we shall assume it equal to 
unity in the following discussion (since we have used £/ = 1 as energy unit). The following 
sections are now dedicated to an analysis of the classical statistical system given by equation 

PD- 



3 Monopoles in dimension 0+1 

We deduce from the interaction (|16D that the interaction betwen two monopoles is 

, . Nt r du r d 2 k cos (k.r) — cos (k.r + lot) . . 

(r ' T) = 12^ J 2^ J (2^p \u)\ (\u\ +7|A;| 3 ) ' ( ' 

We have shifted the interaction by an infinite constant so that the pair interaction between 
two monopoles is finite. The energy of a configuration of monopoles satisfying the neutrality 
condition / q(r, T)dT = is finite and does not depend on the regularization of the pair 
potential. 

From now on, we are interested in the quantum problem at zero temperature, so the system 
is infinite along the imaginary time direction. We shall now use f3 = Nt/12ir to denote the 
inverse fictitious temperature of the monopole plasma, which shouldn't introduce confusion. 
The two parameters associated to ( pT8[) are 7 = l/6k F and the prefactor (3. The interaction 
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( [TE| ) decreases as the distance |r| between two plaquettes increases. In order to have an idea 
of the interaction ranges, we calculate the interaction V(r,r) as a function of r for different 
values of the interplaquette distance. We first rearrange the expression (|i~8|) using the change 
of variables u = ujt and q = (7r) 1 / 3 /c. We obtain 

1 r 1 / 3 / \r\ \ 

with 

. f+°° qdq T . . r+°° du (1 — cosw) 
F x = / ^Jo(gx) / — 1 3 ; (20) 

Jo 27r J-oo 2ir \u\{\u\ + q a ) 

We plotted on figure ^ the interaction for different values of the interplaquette distance \r\. 
In a first approximation, we take into account only the one-plaquette interaction along the 
time direction. In section ||, we renormalize the bidimensional problem with a cut-off for the 
distance between two plaquettes. 

The interactions of the one-plaquette problem are simply V(r) = — r 1//3 F(0)/7 2//3 . We 
look for the phase diagram of the potential V(t) ~ —r a in one dimension as a function of the 
exponent a. Fortunately, some exact results concerning the phase diagram of one-dimensional 



systems with long-range interactions are available [[12]. The main result shows rigorously the 
existence of a finite temperature phase transition for the ID ferromagnetic Ising model if the 
coupling J(n — n') oc l/|n — n'| 7 , with 1 < 7 < 2 []13| . To some extent, these results can be 
transposed to generalized Coulomb gas models, by considering a representation of the Ising 
model in terms of kink and antikink configurations. The potential energy for a single kink- 
antikink pair is then proportional to |n — n'| 2-7 , n and n' being the locations of the kink 
and antikink. The Ising and corresponding Coulomb gas problems are however not equivalent 
since the Ising model generates only rather special configurations where kinks and antikinks 
alternate. We expect intuitively the unrestricted Coulomb gas to be less ordered than the 
corresponding Ising model. By ordered state, we mean the dipolar phase. As a result, an 
unrestricted generalized Coulomb gas with V(r) oc r a is expected to have a high temperature 
Debye-screening phase if a < 1. The fact that a = 1 (the ID genuine Coulomb potential) is 
the borderline is confirmed by several exact investigations [[Uj |15|] showing that this system is 



always in the dipolar phase at any temperature. Our problem is a special case, with a = 1/3. 
We shall now attempt to estimate the transition temperature to the Debye phase. It is then 
tempting to use a real space RG analysis along the lines of references [|l6j-|jl8[. For instance, 



the Coulomb potential in any dimension d (a = 2 — d) has been analyzed in reference [|Tg 
For d > 2, the system is always in a Debye screening phase, whereas for d < 2, there exists a 
finite temperature transition. We note that this simple RG analysis still predicts a non trivial 



fixed point for d = 1 and a = 1, in discrepancy with the exact results of references [O] and 
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T5| . But as d is decreased from 2 to 1, the unstable fixed point is found for higher values of 
the plasma fugacity, so that the dilute approximation leading to the RG equations is no longer 
valid. Hopefully, a = 1/3 is not too large, so the usual RG procedure is consistent. 

In order to analyse the one plaquette problem, we wish to treat the more general problem 
of the generalized Coulomb potential V a in one dimension. We show that if a < 1, the plasma 
has a Debye phase. We call Z T the partition function of the plasma with a minimal separation 
r betwen the charges, which position is allowed to vary from x = to x = L. For a neutral 
system of 2n particles, this defines an integration domain denoted D2 n (L,dr). We wish to 
perform one renormalization step, that is to express Z T as a function of Z T+ s T . To do so, we 
write Z T under the form 



Z T 



E 

m=0 
P 

n 



K 



2m m 



[mi 



E 

p=0 



m 



p\(dr) p / dr 1 ...dr 2 t m -p)W B {ri, 

V) JD 2{m _ p) {L ) T+dr) 



?"2(m-p)) (21) 



(ipi2cosh (PqTE(pi)). 



We have introduced a fugacity denoted by z = Kr. In this equation, we have taken into 
account p dipoles with their center of gravity located at pi (i = l,...,p) and with a size 
between r and r + dr. Wb is the Boltzmann weight of the 2{m — p) remaining isolated charges 
located at ri, r 2 ( m _ p ). E(p) is the electric field at p, created by the 2(m— p) isolated charges. 
More specifically, 



with 



and 



W B (ri,...,r 2n ) = exp ( -(3 1 E < Mi^,ri 



n - r,- 



a ( 


r 






T 





E{ P ) = -Y J qN P v a , T {p 



E 



P-Ti 



a-1 



(22) 



(23) 



(24) 



The p dipoles are assumed to be independent. The integration over the dipole coordinates pt 
must take into account the position of the other charges located at n, ...,r 2 ( m _ p ). We expand 
the cosh in (|2lf) up to second order in the electric field, and we write 



/ 2 cosh ((3qrE(p))dp = 2L + (p(ri, r 2r , 
Jo 



(25) 



where we have used the notation n = m — p and where the integration domain takes into 
account the presence of a hard core condition. We first need to determine the function ip. To 
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do so, we write 



/ 2 cosh ((3qrE(p))dp = E 0(xi+\ — Xi — 3r) 



(26) 



2n 2n 



/ dp 2 + /3 2 rV £ E <&9*K(P - ^)K'(P - * fc ) 

A l+ 2r/2 ^ j=lfc=1 

If we take only the two-body interactions, and the thermodynamic limit, the expression of tp 
takes the form 



f(ri,...,r 2n ) = -3r(2n) + lim /?V 



sign((p - n)(p - Tj))0 



/ rf P 
~T JO 



a-1 



p-r-j 



(27) 



p - n 




p-r-j 


-1) 


T 




T 





+ £? 2 / d P 



p-ri 



T 



2(a-l) 







P - n 



T 



The L — > +oo limit exists provided the system is neutral and a < 3/2. Indeed, the 2n charges 
create a dipolar field at large distances which decays as p a ~ 2 or faster. Taking the square gives 
the upper bound on a for long distance convergency. This property also enables us to shift 
variables and recast the previous expression as a sum of pair contributions which all converge 
separately. We thus obtain 



<p(7"i,-,r 2 „) 



-3r(2n)+p 2 q 2 J2m / d P 

& J -°° 



+oo 



p-ri 



a-l 



p-r 3 



a-l 



(28) 



sign((p - n)(p - rj))B 

2(a-l) 



p-ri 



T 



he 



p-r 3 



P _ n + rj 

r ' 2r 



e 



p 



n + r j 



3 

2 T 



After some rather simple calculations, and extracting the dominant behavior, we have 

2q-1 



<p(ri,...,r 2n ) = f3 2 q 2 T^q i q j c(a) 



Ti - r, 



T 



2nr (3 + (c(a) + d(a))(3 2 q 2 ) . (29) 



The coefficients c(a) and d(a) are given in terms of the Euler B function: 

c(a) = 2 [ a ~ ^ B(a, 2(1 - a)) - B(a, a) 



2a - 1 



d(a) 



2a- 1 V2 



2a-l 



(30) 
(31) 
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We note that (p has the dimension of a length, so that Z T is dimensionless. The scaling 
equations are now obtained and read 



d\np 
dlnr 
dlnK 



a + 2(2a - l)c(a)K 2 r 2 pq 2 



3 + (c(a) + d(a))/3 2 q 4 ) K 2 r 2 



dlnr 2 
And the interaction function is modified by 



dV 
dhir 



2{2a - l)c(a)K 2 T 2 (3q 2 



2a - 1 



2a- 1 



- 1 



A 



T 



1 



(32) 
(33) 

(34) 



where A is the particle separation. These equations are obtained by imposing the normaliza- 
tion constraints V(r) = and 

3S<A-r)~l. (35) 
Since the fundamental form of the interactions is preserved only for the Coulomb potential 
(a = 1), these prescriptions are meaningful mostly near a = 1. In the case a = 1, we recover 
Kosterlitz's RG equations as derived in We note that the second term in the r.h.s. of 
equation fl33|) is not given in |Tj|, but it doesn't change the critical behavior. Its meaning 
is a natural reduction of fugacity because of excluded volume effects. The structure of these 
equations, and in particular, the fact that a > and 2(2a — l)c(a)K 2 r 2 /3q 2 < shows that 
the model keeps a finite temperature transition, and that the size of the Debye screening 
phase increases as a decreases. Coming back to our problem, a is fixed to 1/3 for the single 
plaquette problem. The interaction strength is larger if N increases and if 7 decreases, so 
if kf increases. The dipolar phase is then expected at large N and large electron filling, in 
agreement with the physical intuition that spin-charge separation is more likely to occur in 
the vicinity of the Mott insulator and in the large N limit. This defines a critical Fermi wave 
vector kj(N) such that spin-charge separation occurs for kf > kp(N). The aim of section |4] is 
to obtain quantitative results for the variation of k%(N) with the number of colors N. 



4 Monopoles in dimension 2+1 

We now consider the 2+1 dimensional problem. The derivation of the RG equations is a 
straightforward generalization of what has already been done in the 0+1 dimensional case. 
The potential V(r, r) is given by equation ( jl9|) and (|20|). We use periodic boundary conditions 
in time, with a period L. This regularization will also be used in the numerical calculations. 
Small dipoles (with a length between r and r + dr) can only be parallel to the temporal 
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direction since we require the local neutrality condition / q(r, r)dr = 0, so that the cut-off r 
is only introduced in the temporal direction. The function (p is given by 

pL/2 

(p{rit 1 ,...,r 2n t 2n ) = 2/3 2 g 2 r 2 V V g^fc / dTd o V(0,T)d o V(r k - r 4 - r, A j k - r) 

+2/3 2 g 2 r 2 £ £ g 2 (d V(r, r)f dr. (36) 

r j J3t/2 

In this expression, the summation over r is a summation over the plaquettes which contain 
the small dipoles (with a separation in the temporal direction between r and r + dr). Aj tk is 
the difference between the time coordinates of the monopole j and the monopole k. As in the 
expression fl27|), the integration over the time coordinate of the small dipole contains a hard 
core condition. We evaluated numerically the integrals in (|36|), and took into account only the 
potentials V(r,r) such as |r| < A, with A a lattice cut-off. The RG trajectories are plotted 
on figure (0) for different values of the Fermi wave vector. Notice that some trajectories are 
free to cross each other since the potential depends on the initial conditions via 7. We can 
check the validity of the predictions of the 0+1 dimensional approach: if kf < k^, the plasma 
is deconfined whereas it is confined if kf > k*f. For a one-color model, we find kj = 0.5 ± 0.05. 
However, the Fermi wave vector is bounded above by \/2ti since there is less than 1/2 electron 
of a given color per plaquette. We conclude that there exists a transition even for the one-color 
model. We now adress the question of the N colors monopole model. The action is simply 
multiplied by N, inducing a change in the initial conditions of the renormalization procedure. 
We plotted on figure the critical Fermi wave vector kj as a function of the number of colors. 
We see that for N = 2, which is the case of physical interest, that a possibility of a non Fermi 
liquid arises as the doping increases. 



5 Conclusion 

To conclude, the main result of this investigation is the possibility of tuning the microscopic 
parameters of the model (here the number of colors and the filling factor) in such a way that 
confinement between spinons and antiholons arises or not, depending on these parameters. Us- 
ing a different approach, it had been previously claimed that the Ioffe-Larkin type of plasma 
relevant to the t-J model is always in the dipolar phase, so that spinons and holons have a 
chance not to form a Fermi liquid [10|. It is true that our plasma action (equation (0)) is 



obtained from the simplest fermion loop, without dressing the fermion Green's function, and 
this may be the source of the difference between the results. We were guided here by the very 
strong anisotropy of the intermonopole potential, making it much stronger along the time 
direction, and requiring the neutrality constraint for each plaquette along the time direction. 
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Our intuition is that screening may only be more effective if the two spatial dimensions are 
added to the one plaquette problem, thus weakening the strength of the large distance inter- 
monopole interaction. We hope that these ideas may lead to a more rigourous approach, and 
possibly Monte-Carlo studies of this plasma. Moreover, we have also presented arguments 
showing that a satisfactory microscopic derivation of the plasma action is still missing. One 
difficulty is connected to the ambiguity present while implementing the necessary periodicity 
requirement from a perturbative calculation which by essence assumes 'flat' field configura- 
tions. A second one, and maybe related to the previous remark, is the diverging bare fugacity 
which results from the Villain-type treatment developed in the Appendix. 

B.D. would like to thank J. Wheatley and D. Khveshchenko for their stimulating partici- 
pation in earlier attempts to address the question of collective effects in the monopole plasma, 
and S. Sachdev for an interesting discussion. 



6 Appendix: derivation of the monopole fugacity in the 
Ioffe-Larkin approach 

We start from the quadratic action in the transverse subspace 
1 T 

S = oW 53 H ^) e -L( fc > ^)e±(-/c, -u) + fib{k, u)b(-k, -u). (37) 

In this expression, h is a scalar field and e± is related to b by the Faraday equation 

ub(k, to) = k x e±(k, oj).z. (38) 

Suppose we consider an instanton located on a given plaquette r at time r . The idea is to 
replace b(r, r) by b(r, r) — 27T0(b(r, r) — 7r)5 r ro . Minimizing over transverse configurations of b 
leads to the instanton profile 

G ^ U) = - \ 1 (kZ- 2 + (39) 

b(k, U ) = f-l-H , 2 W)- (40) 

\ to €i(k,uj)uj z + fik z J 

Here, q(k,ou) = 2n exp (— i(k.r + ut )) is the corresponding topological charge density. The 
fugacity is obtained from considering quadratic fluctuations around the single-instanton solu- 
tion and integrating them out. We have to single out the zero mode which corresponds to a 
global translation along the time direction of this solution. This is a standard procedure, and 
we just quote the result fl9|] 

k a(M1X' 2 . (41) 
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In this formula, A is the norm of the zero mode function 



(r,T)~ a -^f±. (42) 



From equation ([40]), we obtain 



It k* V' 2 

A = 2 ^[nZ(^W^W) ' (43) 

The second term is related to the ratio of the product of eigenvalues of the Hessian matrix 
in the vacuum and in the presence of the instanton. In the denominator, it is necessary to 
exclude the zero eigenvalue coming from translation invariance. The Hessian matrix is found 
by expanding the action (|TF|) with the shift b(r, r) — > b(r, r) — 2ii9(b(r, r) — tt), up to quadratic 
order in field deviations 5b(r, r) around the instanton solution. Written in Fourier space, the 
quadratic part of the action is 

62 S = \^Y,^^ h (k^m-k,-uj) (44) 

s k,ui 

-StSEE e^'^-^db^^-k^-u 1 ). 

l ft 0l k ,u k',w' 

We have used the notations A(k,uj) — fi + ei(k,u))u) 2 /k 2 , and 



&o = fWo), (45) 



where b is the instanton profile. The eigenmodes corresponding to equation fl44|) are obtained 
from a rational secular equation since the scattering potential is separable. This equation 
reads 

_ 2 _^Ly i = i. ( 46) 

\b' \N s ^e-A(k,u) 1 ] 

From the structure of this equation, it is possible to derive the determinant ratio as 

Using equation pi]) , and the expression for A leads to K = (fi\b' \Y^ 2 . From equation fl40|), 
we finally get 

— It e V' 2 

k = ^(wX«m*T&) ■ (48) 
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Figure captions 



Figure 1: 

Interaction potential V(|r|,r) of equation (|l~8|) for kf = 4 and for different values of the 
interplaquette spacing as a function of the time coordinate. The notation (0,0) corresponds 
to |r| = 0, (1,0) stands for |r| = 1, (1,1) stands for |r| = v2 and (2,0) stands for |r| = 2. 



Figure 2: 

RG trajectories for N = 1 and different values of the Fermi wave vector. The time coordinate 
is compactified on a circle of length 400. We took r = 1. Some trajectories cross each other. 
This is due to the fact that the potential (|TJ) ( pH|) depends explicitely on 7 and thus on the 
initial conditions. We have plotted the square of the fugacity fugacity z 2 = K 2 t 2 as a function 
of the effective inverse fictitious temperature of the monopole plasma. The dipolar phase 
corresponds to the fixed point (/?, z) = (00, 0) and the Debye screening phase corresponds to 
the fixed point (/3, z) = (0, 00). 



Figure 3: 

Critical Fermi wave vector as a function of the number of colors. The errorbars indicate the 
precision in the location of the fixed point. The dashed line indicates the maximal value of 
the Fermi wave vector (y/2n). 
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